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BOILING  HEAT  TRANSFER 


S.  S.  Kutateladze 


Introduction 


Boiling  heat  transfex-  plays  a  very  important  part  in  the  solu¬ 
tion  of  a  number  of  problems  in  new  technology  and  has  great  appli¬ 
cation  to  industrial  power  engineering.  However,  the  forcing  levels 
of  the  heat-transfer  surfaces  in  these  cases  can  differ  by  several 
orders  of  magnitude.  For  example,  in  refrigeration  engineering, 
when  freon  boils,  the  heat  flirx  has  densities  of  the  order  of 
10^  kcal/m^hr;  in  the  evaporators  of  hydroelectric  stations  they  are 
of  the  order  of  lO^j  in  water -moderated  water-cooled  reactors  of 
nuclear  power  plants  they  are  of  the  order  of  10®;  in  the  elements 
of  Jet  engines  they  are  of  the  order  of  kcal/m^hr. 

In  nuclear  reactors  and  other  highly-forced  Installations  the 
decisive  factor  is  the  critical  density  of  the  heat  flux;  this  deter¬ 
mines  the  upper  safety  limit  for  the  aggregate. 

In  the  apparatus  used  in  industrial  power  engineering  the  abso¬ 
lute  values  of  the  heat-transfer  coefficients  are  of  Importance, 
since  in  such  devices  these  coefficients  are  usually  of  the  same 
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order  for  warming  and  for  heated  media. 

Film  boiling  occurs  In  hardening  processes,  at  near-critical 
pressures,  and  In  certain  other  specific  cases. 

This  report  Is  devoted  to  the  problems  of  heat  transfer  In 
nucleate  boiling.  Nucleate  boiling  is  characterized  by  the  formation 
of  vapor  at  Individual  sites  on  the  heating  surface,  the  so-called 
evaporation  centers,  and  Imparts  to  the  liquid  and  vapor  mixture  a 
very  complex  Inhomogeneous  structure.  This  phenomenon  pertains  to 
one  of  the  most  complex  problems  in  the  hydrodynamics  of  llquld- 
vapor  systems,  namely  to  flows  In  which  the  phase  components  of  the 
flux  are  separated  into  individual  formations  bounded  by  the  mobile 
phase  Interfaces.  These  formations  (bubbles,  drops,  films),  variable 
in  space  and  time,  are  very  numerous,  so  that  here  the  probability 
laws  of  a  multi-element  system  should  be  in  effect. 

At  present  there  are  no  consistent  analytical  methods  for  such 
systems.  The  purpose  of  the  present  work  Is  to  discuss  some  of  the 
main  questions  of  the  boiling  heat-transfer  theory  and  to  generalize 
the  results  of  experimental  v/ork. 

Nomenclature 


C  -specific  heat 

D  =  2R-dlameter 

g  -acceleration  due  to  gravity 

A1  -difference  between  the  mean  enthalpy  of  the  fliuc  and  the 

enthalpy  of  the  saturated  liquid 

L  -length 

p  -pressure 

q  -heat-flux  density 
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q  -critical  heat-flux  density 

^cr 

r  -latent  heat  of  evaporation 

T  -temperature 

T"  -saturation  temperature  above  a  plane 

Wo  -circulation  rate 

X  -vapor  content  of  the  flux  by  weight 

a  -he at- transfer  coefficient 

P  -vapor  content  of  the  flux  by  volume 

7  -specific  weight 

6  -transverse  linear  dimension  of  the  channel 

^  -heating  of  the  flow  core  to  below  the  saturation  temper¬ 

ature 

X  -coefficient  of  thermal  conductivity 

V  -coefficient  of  kinematic  viscosity 

a  -coefficient  of  surface  tension 

T  -time 


Indexes!  "-vapor  phase,  1-phase  interface. 


The  Origin  of  Steam  Bubbles 

The  Intensity  of  nucleate-boiling  heat  transfer  is  due  to  the 
Interaction  of  convection  associated  with  the  directed  motion  of  the 
entire  flow  of  liquid,  and  local  convection,  caused  by  the  origina¬ 
tion,  growth,  and  flotation  of  steam  bubbles.  With  sufficiently 
high  heat-flux  densities,  l.e.,  at  considerable  vapor  formation  rates, 
the  last  process  Is  decisive. 

The  probability  of  the  origination  of  a  vapor-phase  nucleus  In 
a  superheated  liquid  Is  proportional  to  the  number  e_  raised  to  a 
power  which  Is  Inversely  proportional  to  (At)^.  Consequently,  vapor 
bubbles  should  occur  In  the  Immediate  vicinity  of  the  heating  siu’face. 
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where  superheating  of  the  liquid  is  greatest,  compared  with  the  satu¬ 
ration  temperature  above  a  plane. 

Certain  works  on  boiling  heat  transfer  have  used  G,  N,  Kruzhllin’s 
erroneous  viewpoint  [21,  22]  that  vapor  bubbles  originate  on  rough 
projections.  At  present  we  can  state  that  the  boiling  centers  on  a 
solid  surface  are  depressions,  pores,  concavities  at  the  bases  of 
tubercles,  and  rough  places.  Andreyev  [bJj  examining  cylindrical 
pores,  was  first  to  indicate  such  a  mechanism.  This  was  then  examined 
by  iTesis  [33]  and  Courty  and  Foust  [57] »  The  origination  of  a  vapor 
bubble  or  a  conical  pore  was  examined  by  Courty  [57] j  Bankoff  [54], 
and  others,  Zuber  [73]  has  published  a  review  of  these  works. 

Thus,  the  boiling  process  Is  due  not  only  to  the  probability  of 
nucleatlon  with  given  superheating  of  the  liquid,  but  also  to  the 
probability  of  the  distribution  of  vapor-formation  centers  on  the 
heating  surfaces.  The  probable  occxjrrence  of  vapor  bubbles  on  a 
solid  surface  is  clearly  proved  by  the  frequency  distribution  curve 
of  their  formation;  these  were  obtained  for  several  pressures  by 
Zyslna  and  the  author  [15].  Analogous  curves  were  later  obtained  by 
Treshchev  [48]  for  boiling  when  the  liquid  flow  core  Is  heated  to 
the  satiu?ation  temperature.  The  detachment  diameter  of  bubbles  Is 
also  a  variable  value,  subject  to  a  certain  distribution.  Under 
static  conditions,  the  detachment  diameter  of  a  bubble  on  a  smooth 
surface  is  defined,  according  to  calculations  by  Fritz  [60],  by  the 
formula  [23] 

The  variability  of  Dq  Indicates  that  the  value  of  the  edge  angle  0° 
cannot  be  regarded  as  a  definite  characteristic  of  a  great  expanse 
of  the  heating  surface.  In  addition,  recent  motion  pictures  made  at 
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the  Thermophysic B  Institute  show  that  with  Intense  boiling ^  the  bubble 
before  detachment  has  complex  form  and  Is  attached.  In  part,  to  the 
heating  siirface  by  a  cylindrical  neck. 

These  facts  clearly  demonstrate  that  the  hydrodynamic  and  thermal 
processes  during  nucleate  boiling  can  be  described  only  statistically. 


The  Rate  of  Growth  of  a  Vapor  Bubble 


The  rate  of  growth  of  a  vapor  bubble  Is  associated  with  the  heat- 
fl-ux  density  at  Its  surface  by  the  equation 


dR  91 

cfr  ~  %  r"  ‘  ( 2) 

Boshjakowlc  [55]  used  this  relation  to  determine  the  value 
from  the  experimental  data  of  Jakob  et  al.  [64,  65]. 

If  we  consider  the  liquid  medium  around  the  bubble  to  be  Immo¬ 


bile, 


(5) 


where  the  proportionality  factor  C,  according  to  Plesset  and  Zwlck 
[67],  Isy/Jj  while  according  to  Forster  and  Zuber  [59]  it  Is  — ^  l.e., 
C  w  1.6.  Here  T  Is  the  temperature  of  the  superheated  liquid  and 
Is  the  saturation  temperature  at  the  boundary  of  the  bubble. 
According  to  this  formula,  for  the  conditions  of  the  formation 
of  one  of  the  vapor  bubbles  photographed  by  Zyslna  and  me,  we  have 


T,  sec 

0,00475 

0.0095 

0.0285 

q^,  kcal/m^hr 

205,000 

145,000 

85,000 
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Under  conditions  of  thermodynamic  equlllbrlimn,  for  a  spherical 


bubble 

t" -t'LMC- 
'1  '  R(r-rV  cfp 

(^) 


During  heat  transfer,  thermodynamic  equlllbrliim  Is  absent,  the 
bubble  grows  due  to  evaporation  of  the  liquid  In  It,  and  the  valiie 
differs  from  ( 4) ,  Assuming  tha^  there  Is  only  vapor  In  the  bubble 
we  can.  In  the  first  approximation,  relate  the  evaporation  rate  and 
the  pressure  drop  In  the  vapor  boundary  layer  by  the  familiar  kinetics 
formula 


(5) 


where  ^  <  1  Is  the  accomodation  coefficient,  B  Is  the  universal  gas 
constant,  and  m  Is  the  molecular  weight. 

The  corresponding  change  in  the  saturation  temperature 


(6) 

Adding  ( 4)  and  ^  6)  .  we  find  that 


where  ^  In  technical  units,  A  =  4^^"* 

Figure  1  compares  the  terms  In  brackets  in  Formula  (7)  •  From 
the  figure  It  is  clearly  evident  that  even  when  ^  =  1  the  correction 
associated  with  evaporation  kinetics  Is  not  only  commensurate  with 
the  correction  for  the  curvature  of  the  phase  Interface,  but  in  a 
number  of  cases  can  even  exceed  It. 
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Pig.  1.  Values  of  AT"  and  AT"  according 
to  Pcrmula  (7). 

First  Heat  fflux  Critical  Penalty  with  Flow  In  Tubes  and  Channels 

Any  examination  of  the  boiling  process  must  begin  with  the 
problem  of  free  convection  In  a  large  volume  of  liquid.  However,  In 
technology  the  most  essential  problem  is  the  forced  flow  of  a  boiling 
liquid  in  tubes  and  channels. 

At  the  present  time  much  experimental  data  are  available  on 
critical  heat  fluxes  in  tubes,  while  there  Is  not  much  available  on 
fluxes  In  slit  channels. 

It  is  extremely  important  to  explain  the  effect  of  flow-rate 
pulsations  on  the  value  q^j^.  There  are  a  great  many  operational  and 
therraophysical  parameters  which  determine  the  limiting  density  of 
the  heat  flux  for  which  nucleate  boiling  and,  consequently,  high 
heat-transfer  intensity  are  still  maintained.  In  this  regard,  the 


question  of  methods  of  a  generalized  analysis  of  the  examined  phenom¬ 
enon  and  rational  forms  of  the  connection  between  its  determining 
values  is  of  priraiary  Importance. 

Two  hypotheses  have  been  advanced  with  regard  to  the  boiling 
crisis  in  large  volumes  of  liquids. 

The  first  of  these^,  proposed  by  the  author  [24^  25]  j  presumes 
that  the  transition  from  nucleate  to  film  boiling  is  a  qualitatively 
specific  process j  the  basis  of  which  is  the  loss  of  stability  of  the 
structure  of  a  two-phase  boundary  layer  which  existed  prior  to  the 
crisis , 

In  this  case,  as  for  any  hydrodynamic  crisis,  the  transition 
from  one  boiling  point  to  another  is  determined  not  by  a  strictly 
fixed  value  (in  this  case  q_„)  ,  but  only  by  its  most  probable  value. 

V'X 

The  second  hypothesis,  most  clearly  formulated  by  G,  N,  Kruzhilin 
[21]  is  based  on  the  assumption  that  during  the  first  boiling  crisis 
no  new  qualitative  peculiarities  arise  and,  consequently,  the  critical 
heat  flux  density  is  determined  by  the  same  initial  group  of  general¬ 
ized  parameters  (  simllairlty  criteria)  as  is  the  heat-transfer  coeffi¬ 
cient  during  nucleate  boiling.  This  hypothesis  has  not  been  developed 
for  forced  flow,  while  for  free  convection  boiling,  as  has  been 
pointed  out  [28],  it  has  resulted  in  an  unjustifiably  complex  empir¬ 
ical  formula. 

According  to  the  hydrodynamic  hypothesis,  the  transition  from 
nucleate  to  film  boiling  begins  with  a  less  in  the  stability  of  the 
two-phase  boundary  layer  in  one  of  the  small  regions  of  the  heating 
surface,  and  the  formation  of  a  stable  vapor  "spot."  The  increase 
in  temperature  occurring  beneath  this  "spot,"  due  to  a  decrease  in 
the  heat-transfer  coefficient,  results  in  further  development  and 


-8- 


propagation  ol  film  boiling  along  the  heating  surface.  Such  a  system 
is  in  complete  agreement  with  that  observed  during  boiling  on  heat- 
resistant  (e.g.j  graphite)  rods  and  plates. 

It  Is  not  too  difficult  to  describe  the  conditions  of  the  mechan¬ 
ical  interaction  between  vapor  and  a  liquid  for  a  differential  element 
of  the  phase  Interface  of  such  a  system.  Here,,  since  it  is  a  ques¬ 
tion  of  the  local  stability  of  structure,  the  matter  of  boundary 
conditions  for  the  entire  system  does  not  arise.  Accordingly,  the 
system  of  equations  consists  of  motion  and  continuity  for  each  phase, 
and  conditions  of  the  mechanical  phase  Interaction  at  the  Interface 
[25]. 

Treating  an  Invlscld  liquid,  we  can  also  use  Rayleigh's  theory 
of  the  stability  of  a  double-layer  flux  [39], 

When  the  problem  is  treated  thus,  the  value  to  be  determined  is 
the  critical  vapor  rate,  and  the  question  arises  as  to  the  connection 
between  the  local  value  and  the  mean  heat  flux  density  at  the  heating 
surface . 


Evidently,  the  only  actual  measure  of  the  rate  of  vapor  formation 


Is  the  value 


the  maximum  amount  of  vapor  which  can  form  per  unit  of  the  heating 
surface  per  unit  time. 


Here  we  will  examine  the  maximum  amount  of  vapor,  since  in  the 
general  case  part  of  the  heat  flux  q  may  not  be  expended  on  vapor 
formation  (e.g.,  In  a  flux  with  a  core  which  is  heated  to  below  the 
saturation  temperature) , 

Yet  another  characteristic  Is  evident  from  the  heat  balance  of 


the  mixture;  this  can  be  represented  as  the  criterion 
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where  A1  is  the  difference  between  the  mean  enthalpy  of  the  flux  and 
the  enthalpy  of  the  saturated  liquid. 

1 

In  a  flow  of  satirrated  liquid,  =  x,  where  x  Is  the  vapor 
content  by  weight,  l.e..  In  this  case  the  criterion  K  actually  has 
purely  hydrodynamic  significance. 

When  ![.'  <  T",  the  criterion  K  =  ■  where  =  T"  -  T  Is  the 

subheating  of  the  flow  core  to  the  saturation  temperature.  In  this 
case,  K  Is  purely  hydrodynamic  only  when  the  heat-transfer  mechanism 
Is  completely  turbulent]  this  Is  a  measure  of  the  decrease  In  vapor 
flux  due  to  condensation  In  the  cold  core.  If  the  molecular  thermal 
conductivity  Is  significant,  the  criterion  K  Is  associated  with  the 
Prandtl  number  Pr  as  a  thermal  criterion.  We  will  not  discuss  this 
at  this  time. 

The  stability  of  a  two-phase  system  is  always  influenced  by  the 
kinetic  head  of  the  originating  phase,  the  measure  of  which.  In  this 
case.  Is  the  value 

il' 

i  ’  (10) 

gravity,  and  surface  tension.  As  the  measure  of  the  Joint  action  of 
these  two  latter  forces  ( keeping  (1)  in  mind)  we  can  use  the  value 


(11) 

The  square  root  or  the  ratio  of  these  values  Is  the  criterion 
[24,  25] 


When  q  =  q^^  criterion  (12)  becomes  definite. 


(12) 


If  we  use  the  Rayleigh  stability  criterion,  criterion  K  assumes 


the  form 


(15) 

The  value  of  this  group  for  the  first  critical  heat-flux  density 
(the  occurrence  of  film  boiling)  is  designated  by  the  Index  1;  for 
the  second  critical  heat-flux  density  (  destruction  of  the  vapor  film) 
the  value  Is  designated  by  the  Index  2, 

The  difference  between  expressions  ( 12)  and  ( 13) >  the  value 
"sjl  j  has  no  essential  significance,  and  Is  associated  with  the 
solution  of  the  Initial  system  of  equations  used  by  Rayleigh. 

Various  authors  [24,  25^  26]  have  shown  that  within  the  frame¬ 
work  of  these  assumptions  we  have  the  following  connection  between 
the  dimensionless  criteria; 

k yr J  >  ^  I  ^  y  ’  r'l  J  T'  'tj  (14) 

Here  we  should  give  a  distribution  law  for  the  feed  of  heat 
along  the  heating  surface. 

For  an  underheated  liquid  the  hydrodynamic  hypothei^ls  has  re¬ 
sulted  In  a  conclusion  as  to  the  existence  (  at  least  In  first  approx¬ 
imation)  of  a  linear  relationship  between  the  criteria  K  and 
In  the  form 


q 


cr  rf 


(15) 


For  purely  turbulent  transfer  the  coefficient  B  Is  a  certain 
constant.  If  molecular  heat  conductivity  Influences  the  heat  trans¬ 
fer  near  the  nascent  elements  of  the  vapor  phase  [e,g,,  according 
to  (  3) ] ,  the  coefficient  B  Is  a  function  of  the  Ar  and  Pr  criteria. 

Turning  to  experimental  data,  let  us  note  that  despite  numerous 
Investigations j  the  various  results  are  not  In  satisfactory  agree- 
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merit.  Here  we  are  speaking  not  only  of  the  purely  quantitative  dis¬ 
crepancies  which  can  be  explained  by  the  instability  of  the  data  on 
critical  regimes. 


A  number  of  the  discrepancies  were  explained  after  the  role  of 
preincluded  volumes  in  the  occurence  of  pulsation  and  the  Influence 
of  the  length  of  the  heated  segment  had  been  explained  [4,  11 ,  46], 
Most  experiments  on  slit  channels  and  tubes  indicate  that  the 
influence  of  the  slit-width  or  the  diameter  is  slight ^  when 


The  linear  relationship  between  q  and  Ai  in  the  underheating 
region  >  0  is  satisfactorily  verified  by  experimental  data  for 

relatively  open  channels  with  pressures  up  to  critical,  Povarln  and 
Semenov  have  verified  the  satisfactory  nature  of  this  relationship 
for  liquid  flow-rates  to  4o  m/sec  [36], 

The  experiments  by  various  authors  [12,  13 j  17 ^  ^6,  52]  have 
distinctly  verified  the  general  natixre  of  relation  (15).  However, 
the  numerical  values  of  B  and  n  are  not  the  same  in  various  experi¬ 
ments  , 


The  experiments  of  Zenkevich  [13]  and  Mlropol’skly  and  Shltsman 
[17]  give  a  clear  relationship  between  q^^  and  Wq  when  ^  >  Oj  this 
is  satisfactorily  described  by  the  empirical  formula  [24,  25] 


n  ^0,  its- 


(17) 


where  Kio  is  the  value  of  the  complex  Ki  when  -S-  =  0, 

The  experiments  of  Chlrkin  and  Yukln,  Averin  and  Kruzhilin, 
and  the  author,  on  the  flow  of  water  in  wide  annular  channels  are 
also  satisfactorily  described  by  a  formula  of  the  type  in  ( I7)  but 
with  a  proportionality  factor  of  0,085  [25,  26],  Ornatskly's 
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results  [35]  are,  quantitatively,  about  50^  lower. 

Experiments  by  Doroshchykln  and  Frid  [11,  1?]  have  showri  the 
weaker  effect  of  liquid  flow-rate  on  q„„,  while  there  is  a  slight 

Ox 

Influence  of  Wq  on  the  coefficient  B, 

The  experiments  by  Alad'yev  and  Dadonov  [3]  differ  in  this 
respect  from  the  data  examined  above.  In  these  experiments  Wq  has 
practically  no  effect  on  q,  „  at  low  pressures,  while  at  high  pres- 

OX 

sures  this  effect  is  relatively  slight  and  the  monotonlc  relation 

y  n 

between  the  coefficient B( -r^i)  and  press'ure  is  destroyed. 

The  work  by  Povarnln  and  Semenov  [36]  contains  a  discussion  of 
the  connection  between  q  and  Wq  for  high  liquid  flow-rates, 

OX 

In  the  region  x>  0,  i,e,,  in  the  flow  of  a  saturated  vapor- 
liquid  mixture,  the  relationships  for  q^  are  even  more  complex.  It 
has  been  reliably  established  that  the  degree  of  Influence  of  vapor 
content  depends  on  the  mass  flow-rate  of  the  mixture  or,  in  other 
words,  on  the  so-called  circulation  rate 

^_aaL 


(18) 


where  is  the  mass  flow-rate  of  the  mixture  per  second,  and  0,  is 
the  cross  section  area  of  the  canal. 

In  the  absence  of  pulsations  l,e.,  the  critical  heat 

flux  density  decreases,  despite  the  essential  Increase  in  the  mean 
flow-rate  of  the  liquid  with  increased  vapor  content  of  the  mixture. 
Comparing  this  result  with  ( IJ) ,  we  must  conclude  that  increased 
vapor  content  of  the  flux  essentially  changes  the  vapor  content  of 
the  near-wall  layer  and  facilitates  destruction  of  the  liquid 
structures  existing  in  this  layer. 

From  the  available  experimental  data  we  can  note  a  more  or  less 
clear  dependence  only  for  the  region  of  vapor  content  by  volme 
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p  <  0,9.  Since  the  vapor  contents  by  volxune  and  by  weight  are  asso 
ciated  by  the  formula 


f  i-X 


(19) 


the  values  p  <  0.9  at  low  pressures  correspond  to  very  slight  vapor 
content  by  weighty  while  at  high  pressures  the  vapor  content  by 
weight  is  of  the  order  of  0,1  and  greater, 

R,  A.  Rybln  has  processed  the  experimental  data  of  Alad’yev, 
Bukhberg,  Doroshchuk,  Dzhekeg,  Layne,  Miropol • skly,  Subbotin,  et  al. 
pertaining  to  tubes  of  9Mm  and  more  diameter ^  the  vapor  content 
0<p<0,9^  and  inlet  conditions  where  the  water  entering  the  exper 
imental  sector  was  heated  to  just  below  the  saturation  temperature. 
The  results  for  two  values  of  the  criterion 


'^0 


(20) 


are  shown  in  Pig,  2,  from  which  can  be  seen  the  absolute  discrep¬ 
ancies  among  various  observations,  and  the  over-all  scatter  of  the 
experimental  points. 

It  is  interesting  to  note  that  the  average  lines  at  point 
P  =  0  practically  coincide  with  calculations  by  Formula  ( I7)  . 

Figure  3  shows  averaged  straight  lines  for  var'lous  values  of 
criterion  (20).  As  can  be  seen,  as  this  criterion  decreases,  i.e,, 
the  circulation  rate  decreases,  the  relative  vapor  content  also 
decreases . 

For  the  region  of  high  vapor  content  by  weight  (P  >  0,9)  there 
are  not  enough  available  experimental  data  for  generalizations. 

Thus,  systematic  compilation  and  further  accumulation  of  exper¬ 
imental  data  Is  a  necessary  condition  for  constructing  reliable  com¬ 
putational  dependences  for  with  flow  in  tubes  and  channels. 


-14- 


r  ■  » 

— ; — 1 — ! — 

TTl 

^-rrr- 

-  .  1  ■  <  1  '  '  1 

xtt 

._U-4 

j  1 

Id 

.  !  }  : 

;  ;  i  ■  : 

•  1' '  1  *  :  : 

-  !i4 - 0— 

: - 

— Lj 

.f 

-  ■ 

Vt  1 
r-  = 

- h-'— V 

-1  •  ! 

Fig,  2,  The  value  Kj  vs.  p  and  P  for  flow 
In  tubes,  p  <  0.9, 


Fig,  5.  The  value  Ki  vs.  p  and  F*  according 
to  averaged  lines. 


In  this  regard,  the  complex  empirical  formulas  proposed  by  Zenkevich 
and  Ivashkevich  [15,  16]  cein  be  considered  as  a  first,  and  not  very 
successful,  attempt. 

Critical  Heat-Flux  Density  with  Free  Convection  In  a 
Large “Volume  Liquid 

The  results  of  applying  the  hydrodynamic  theory  of  boiling 
crises  to  conditions  of  free  convection  In  a  large-volume  liquid  are 
qulle  familiar,  and  reduce  to  the  following  basic  relationships 
[24,  251: 


Kio  =  constj 


(21) 
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<  1. 


(22) 


2 

^cr,  1 


Formula  (  13)  expresses  the  effect  of  not  heating  the  llquild  mass  to 
the  saturation  temperature. 

Using  the  Rayleigh  theory  of  the  stability  of  a  two-layer  flux, 
Zuber  and  Trlbus  [72j  found  that  In  an  Invlscld  liquid 


<  0,1$ 

3 


(23) 


Relation  (  3)  was  also  used  In  calculating  the  latter  two  coeffi¬ 
cients. 

It  Is  Interesting  to  note  that  complex  B,  according  to  (23)  »  Is 
a  very  slowly  changing  value.  The  theoretical  values  of  (23)  closely 
agree  with  most  experimental  data. 

However,  It  Is  Important  not  to  be  concerned  with  this  good 
qualitative  and  quantitative  verification  of  the  theory,  but  to 
discuss  the  various  discrepancies  noted.  Evidently,  If  these  devia¬ 
tions  are  not  associated  with  inaccurate  experiments  and  data  on  the 
physical  properties  of  the  Investigated  substance,  they  should  serve 
as  the  starting  point  for  further  development  of  the  theory.  We  are 
speaking  of  small  (<  0.1)  values  of  Ki  for  molten  salts  using  approx¬ 
imately  the  same  values  for  Yi^  as  in  Robin's  experiments  [38],  and 
deviations  of  B  and  n  In  experiments  for  biphenyl  from  the  values  of 
these  factors  for  water  and  alcohol  [27]» 

Very  Interesting  results  on  critical  heat  fluxes  in  mixtures 
have  been  obtained  by  van  Wljk,  van  Stralen  [70],  and  Fastovskly 
and  Artym  [51]  for  thln-wlre  boiling.  They  found  that  when  a  binary 
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considerably 


mixture  has  a  specific  composition,  the  value  q 
exceeds  the  corresponding  values  for  each  of  the  components.  For 
pure  components  the  obtained  values  for  the  critical  heat  flux  are 
In  good  agreement  with  Formula  (21),  when  Borlshanskly ' s  corrections 
for  viscosity  [6,  25]  are  Introduced; 

■  (24) 

The  mechanism  of  this  phenomenon  Is  far  from  clear.  However, 
the  experiments  conducted  by  Bobrovlch  et  al,,  at  the  Therraophyslcs 
Institute,  using  heating  surfaces  larger  than  those  In  the  above- 
mentioned  experiments,  did  not  reveal  this  effect  In  any  clearly 
recognizable  form  (Fig,  4)  . 


kcal/m^hr 


Fig.  4.  The  value  ^  vs,  concentration 

of  aqueous  butanol;  •  -  wire,  p  =  9*8  ^^tm 
(abs.);  O  -  plate,  p  =  11  atm  (abs.)  . 


The  Effect  of  Forced  Convection  on  Heat  Transfer  with 
Nucleate  Bolling  In  Tubes 

The  Influence  of  directed  convection, characterized  by  the  mass 
flow  rate  and  the  vapor  content  of  the  flux,  on  the  intensity  of 
nucleate-boiling  heat  transfer  depends  on  the  heat-flux  density  and 
the  physical  properties  of  the  heat-transfer  agent,  Sterman, 
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Styushln,  Sorokin,  et  al.,  [18,  hi  -  44]  have  devoted  much  effort  to 
this  problem. 

With  high  pressiice  and  strong  heat  fluxes,  the  flow  rate  depends 
only  slightly  on  the  nucleate-boiling  heat-transfer  coefficient;  this 
has  been  verified  by  the  experiments  of  Alad'yev,  Armand,  et  al,, 

[2,  47], 


At  low  pressure  and  moderate  heat  fluxes  this  effect  can  be 
quite  essential. 

As  the  author  has  shown  [24,  29],  the  effect  of  the  Influence 
of  circulation  rate  on  the  nucleate  heat-transfer  coefficient  can  be 
eatlmated  fr'om  the  ratio  of  the  heat-transfer  coefficients  during  the 
flow  of  a  nonboiling  liquid  ao  and  during  developed  boiling  aoo* 

In  that  region  where  vapor  formation  has  no  effect,  we  have 


—  r  f  I 


(25) 


whereupon 


(26) 


The  simplest  Interpolation  formula  which  satisfies  conditions 


(26)  Is 


(27) 


Figure  5  shows  a  vs.  Wq  and  q  from  one  of  Sterman's  experiments,. 
From  the  figure  It  Is  clearly  evident  that  the  influence  of  circula¬ 
tion  rate  decreases  with  Increasing  £, 

Processing  of  these  experiments  has  shown  that  a  value  n  =  2 
can  be  considered  convenient  for  practical  calculations  [2:9J. 

The  value  aooj  the  heat-transfer  coefficient  during  developed 
boiling,  is  determined  for  a  given  liquid  by  the  formula 
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Pig,  5.  The  value  -5L-  vs.  q  and  Wq 

qO*T 

for  water,  p  =  2  atm  (abs.)  ,  tube 
diameter  l6  mm. 


(28) 


The  values  of  the  pressure  function  C(p)  according  to  experi¬ 
ments  with  the  boiling  of  water  In  tubes  are  given  In  Fig.  6  from 
data  of  Alad'yev,  Sterman,  et  al.  [2,  44,  47],  This  relation  Is 
qualitatively  analogous  to,  but  somewhat  lower  than,  that  obtained 
by  Borlshanskly  [7]  in  his  experiments  with  free-convectlon  boiling 
in  a  large-volume  liquid.  This  20-30^  quantitative  discrepancy  is 
of  relatively  slight  significance. 

This  dependence  holds  when  p  <  0,7.  It  should  be  noted  that  the 
available  experimental  data  are  by  no  means  ample  for  these  vapor 
contents.  For  the  region  p  >  0,7>  however,  the  available  data  are 
insufficient  for  generalizations. 

Experiments  by  Tarasova,  Armand,  and  Kon'kov  [47] ^  conducted  at 
a  pressure  of  I70  atm  (abs.)  have  shown  that  when  p  >  0.7  there 
occurs  first  an  increase  and  then  a  sharp  drop  in  the  value  of  the 
heat-transfer  coefficient. 
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atm  (  abs.) 

ilg,  6,  The  value  C(p)  In  Formula  (28)  from 
experiments  with  boiling  water  In  tubes:  •  - 
according  to  [2]:  O  -  according  to  A  - 

according  to  [^fjJ  curve  1  according  to 
Kruzhllln's  formula;  curve  2  according  to 
Kutateladze 's  formula.  In  the  general  form¬ 
ulas  the  proportionality  factors  are  reduced 
by  20^. 


Developed  Nucleate-Boiling  Heat  Transfer 

The  lack  of  Influence  of  velocity  and  channel  dimensions  on  the 
heat-transfer  coefficient  during  developed  nucleate  boiling  ( at  least 
when  3  <  0.7  and  D  >  5  'f^'  )  is  associated  with  the  fact  that  In 
this  case  the  washing  of  the  heating  surface  by  the  liquid  Is  com¬ 
pletely  determined  strictly  by  the  process  of  vapor  formation.  This 
Is  yet  another  verification  of  the  qualitative  difference  in  heat- 
transfer  mechanisms  during  nucleate  boiling  and  during  the  origination 
of  film  boiling,  since  the  value  q^j,  is  quite  clearly  a  ftmctlon  of 
velocity  in  those  cases  where  a  is  practically  Independent  of  it. 

Thus,  during  developed  boiling,  only  the  heat-flux  density  q 
is  an  operational  factor  with  respect  to  the  heat-transfer  coeffi¬ 
cient,  In  addition,  influence  is  exerted  by  the  group  of  physical 
properties  of  the  medium  v/hich  change  along  the  satiuratlon  line. 


-20- 


The  probable  distribution  of  vapor  formation  centers  along  the 
heating  surface  and  the  process  of  the  generation  of  vapor  bubbles 
cause  the  hydrodynamic  and  thermal  situations  to  be  unstable  In  the 
near-wall  region^  as  has  been  mentioned  previously. 

In  this  regard,  the  values  of  the  heat -transfer  coefficients 
are  far  from  identical  in  various  experiments.  Proper  attention  is 
not  devoted  to  this  clrcurastance ,  and  Individual  experimenters  per¬ 
sistently  consider  that  only  their  data  are  reliable. 

When  determining  the  criteria  of  similarity  of  the  heat-transfer 
process  It  Is  necessary  to  examin''  the  interaction  of  the  entire 
flux  with  the  heating  surface,  l.e.,  the  Investigation  must  not  be 
restricted  merely  to  local  conditions,  as  was  done  in  the  problem  of 
the  critical  regime.  This  causes  as  yet  Insurraoimtable  difficulties 
in  the  forimulatlon  of  a  corresponding  system  of  equations  to  which, 
in  addition,  must  be  added.  In  the  given  case,  conditions  at  the 
boundaries  of  the  system. 

However,  these  difficulties  differ  little  from  similar  diffi¬ 
culties  in  the  analytical  formulation  of  ordinary  turbulent  flow 
where  the  flux  also  consists  of  the  statistical  aggregate  of  final 
formations  which  are  variable  in  time  and  space. 

Since  during  boiling  various  areas  in  space  have  differing 
physical  properties  ( liquid  and  vapor)  ,  the  mathematical  formulation 
of  the  problem  should  encompass  relationships  for  the  phenomena 
within  each  region  and  also  for  the  Interactions  (mechanical  and 
thermal)  on  their  surfaces,  l.e.,  at  the  phase  interfaces* 

The  conditions  on  the  boimdary  of  the  entire  system,  l.e.,  on 
the  heating  surface,  should  be  given  Independently, 

All  reactions  in  a  two-phase  flux  of  any  complexity  are  a  direct 
result  of  a  group  of  local  reactions.  Therefore,  If  we  assume  that 
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this  aggregate  does  not  Introduce  any  particiilar  physical  magnitude 
the  apparatus  for  deriving  the  slrillarlty  criteria  consists  of  a 
system  of  differential  motion  and  thermal-conductivity  equations  In 
each  phase,  the  conditions  of  local  reactions  at  the  Interfaces,  and 
the  conditions  on  the  heating  surface  (proper  boindary  conditions). 

This  problem  has  been  presented  In  such  form  by  the  author 
[23,  24]  and  later  by  Eygenson  [53]* 

This  approach  to  the  derivation  of  similitude  conditions  can 
be  formulated  In  two  ways: 

1)  we  write  tne  heat-transfer  and  motion  equations  In  each  phase 
In  differential  form,  the  equations  of  the  reaction  on  a  differential 
element  of  the  Interface,  and  the  conditions  at  the  system  boindarlesj 

2)  we  write  the  heat-transfer  and  motion  equations  In  differ¬ 
ential  form  for  the  dominant  phase  In  the  process  ( for  the  liquid  In 
nucleate  boiling) ,  The  second  phase  Is  treated  as  the  final  elements 
of  specific  form  (e.g,,  spherical  vapor  bubbles)  ,  We  wi’lte  the 
motion  equations  for  the  center  of  gravity  of  such  an  element  and  the 
heat-balance  equations.  The  conditions  at  the  bouadarles  of  the 
system  are  given  as  In  the  first  case . 

We  can  also  compile  a  system  of  equations  by  considering  the  new 
phase  as  a  heat  and  mass  sink  [10,  42].  Methodologically,  this  Is 
no  different  from  that  given  above,  but  In  this  case  the  discrete 
elements  of  the  new  phase  must  be  considered  as  material  points. 
Therefore,  such  a  method  can  be  used  when  the  region  encompassed  by 
the  process  Is  very  large  compared  with  the  elements  of  the  dispersed 
phase . 

Both  methods  are  basically  identical,  since  the  equation  of 
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motion  of  the  center  of  gravity  of  a  bubble  can  be  obtained  by  solv- 
ing  the  combination  of  motion-equation  systems  for  both  phases. 

Such  solutions  have  been  in  final  form  for  individual  regions  of 
Reynolds  numbers  by  Hadamard  [61],  Rybchlnsky  [69],  and  Levlch  [50]. 

Naturally,  the  system  of  equations  should  correspond  to  the  model 
used  for  the  actual  physical  process,  and  be  closed,  in  the  sense 
that  the  number  of  equations  corresponds  to  the  nimiber  of  unlmowns 
contained  therein.  However,  the  degree  of  correapondence  of  the 
model  to  the  actual  process  is  determined,  in  final  analysis,  by 
comparison  with  the  results  of  generalization  of  the  experimental 
data.  This  factor,  naturally,  is  not  a  feature  of  the  boiling  pro¬ 
cess,  but  has  general  significance  [19]* 

Examination  of  the  motion  and  thermal-conductivity  equations 
of  the  liquid,  taking  into  account  the  equations  of  motion  and  of 
mechanical  and  thermal  reactions  of  the  vapor  phase,  yields  three 
criteria: 

P. _ 1-/  (T 

ir/'V  ‘  ir-ir"}  * ' 

(29) 

Here,  as  the  velocity  scale,  we  have  the  value  ,  while  as  the 

linear-dimension  scale  we  have  Wyy^- 

With  very  Intense  boiling,  the  near-wall  layer  can  be  so  agi¬ 
tated  that  the  Influence  of  molecular  viscosity  can  be  disregarded 
in  the  motion  equation.  In  such  a  case.  Instead  of  the  three  criteria 
(29}  we  have  only  two: 

(30) 


^  rrKe^  • 
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This  holds j  In  particular j  for  media  with  Pr  «  1.  From  the  boundary 
conditions  it  follows  that 


^  (r-W  ■ 


(31) 


This  criterion  was  introduced  by  Jakob  [65], 

The  condition  of  continuity  of  the  temperature  field  at  the 
phase  interfaces  has  the  form 


h,  (32) 

If,  as  was  done  by  Kruzhllln.  [21],  we  define  T”  according  to 
Formula  (4)  we  get  another  new  criterion  which  can  be  represented 
In  the  form 


(vr"/ 

Acrjycir-f") 


(33) 


If  we  consider  the  kinetics  of  evaporation,  as  was  done  by 
Zysln  [14,  24].,  we  have  the  criterion 


As  was  shovm  above. 


(y*) 


the  terras  of  the  equation  that  result  in 


these  two  criteria  are  at  least  conjugate, 

dp  ^  p 

If  we  assume  that  we  have 


(35) 


where 


(36) 


It  is  difficult  to  process  experimental  data  with  the  inclu¬ 
sion  of  the  criteria  and  (or  and  ,  It  appears, 
however,  that  for  many  liquids  there  is  practically  a  single-value 
relationship  between  these  criteria.  This  explains  the  fact  that 
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empirical  formulas  constructed  for  a  limited  system  of  criteria  of 
the  type 

*■,;  )  (37) 

give  more  or  less  satisfactory  results. 

For  a  specific  liquid,  one  of  these  processing  methods  has 

slight  preference,  Tohllevlch  [49]  and  certain  other  authors  have 

shown  that  processing  using  the  criterion  K  gives  better  results, 

]P 

particularly  In  the  region  of  high  pressures.  In  particular,  this 
is  true  of  water,  an  important  working  substance  (Fig.  6). 

In  striving  to  Improve  the  agreement  with  experimental  data, 
certain  authors  have  constructed  formulas  that  contain  a  limited 
number  of  criteria  and  have  Individual  constants  for  each  liquid. 

In  certain  formulas  these  constants  are  even  dimensional.  Such 
formulas  have  been  derived  by  Rohsenow  [68],  Alad'yev  [2],  Jakob 
[66],  Tolublnskly  [50],  and  Levy  [62],  Evidently,  In  this  case 
simple  empirical  formulas  for  a  given  mass  are  preferable, 

Alad'yev  has  given  a  successful  formula  of  this  type  for  water; 
this  can  be  given  conveniently  In  the  form 

'  (  38) 

where  C  Is  a  constant  which  characterizes  the  heating  surface.  With 
tube  boiling,  C  2,1;  for  large -volmie  boiling  In  brass  pipes 
C.;^2.6;  for  large-volume  boiling  In  stainless  steel  pipes  G:3r3. 

Methods  of  constructing  empirical  formulas  on  the  basis  of  the 
theory  of  corresponding  states  of  matter  are  also  useful  In  this 
respect.  This  question  has  been  treated  by  Novikov  [34],  Rychkov 
[40],  Borlshanskiy  [8],  and  others. 
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Conclusions 


1,  The  available  theo'retlcal  and  experimental  data  on  the  mech¬ 
anism  of  nucleate  boiling  pertain  mainly  to  conditions  in  which  the 
bubbles  still  react  with  one  another  only  weakly.  Mien  a  saturated 
liquid  boils  in  a  region  of  high  heat-flux  densities,,  the  mechanism 
of  the  development  of  vapor  bubbles  changes  substantially.  Further 
investigations  of  such  types  are  necessary  to  construct  a  relatively 
complete  theory  of  nucleate  boiling. 

2.  The  hypothesis  of  the  hydrodynamic  nature  of  crises  in 
boiling  and  the  development  of  a  crisis  from  local  disturbances  in 
the  stability  of  the  two-phase  boxmdary  layer  has  made  it  possible 
to  formulate  a  system  of  differential  equations  to  describe  this 
process,  to  select  a  rational  form,  of  the  specific  criterion,  and 
select  a  system  of  basic  specific  criteria. 

On  the  basis  of  this  hypothesis,  theoretical  formulas  have  been 
derived  for  the  first  critical  heat-flux  density  in  a  large  volume 
of  free-convection  liquid. 

The  qualitative  conclusions  from  the  hydrodynamic  hypothesis 
are  also  verified  by  data  on  forced  flow. 

Further  development  of  this  problem  requires  vast  study  of  the 
critical  phenomena  for  the  flow  of  liquids  with  various  physical 
properties  and  accurate  notation  of  the  experiment  conditions. 

It  is  particularly  important  to  gather  experimental  data  on  the 
positive  vapor  contents  of  a  flux, 

5.  The  influence  of  circulation  rate  on  the  nucleate-boiling 
heat-transfer  coefficient  depends  on  the  relation  of  Intensities  of 
heat  transfer  during  ordinary  convection  and  during  developed 
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boiling.  In  certain  cases  this  Influence  can  be  quite  substantial 
and  aid  in  intensifying  the  operation  of  the  heat-exchanger, 

4.  Strict  formulation  of  the  problem  of  the  intensity  of  nucleate- 
boiling  heat-transfer  is  possible  only  by  examining  statistically 

the  multitude  of  final  elements  which  are  variable  in  space  and  time. 
Such  statistics  are  lacking  at  the  present  time . 

The  above-mentioned  methods  for  circumventing  this  difficulty 
in  deriving  similitude  conditions  result  in  a  relatively  complete 
system  of  hydrodynamic  and  thermal  criteria.  Conditions  on  the  heat¬ 
ing  surface  cannot  be  strictly  formulated, 

5.  Practically  no  studies  have  been  made  of  conditions  of  heat- 
transfer  during  the  boiling  of  binary  mixtures  and  for  certain  other 
special  cases.  The  amount  of  work  done  along  this  line  in  no  way 
corresponds  to  the  practical  and  theoretical  significance  of  the 
problem. 
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